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Abstract 

We have performed an ab initio study of the thermodynamical properties of rare-earth-magnesium inter- 
metallic compounds MgRE (RE=Y, Dy, Pr, Tb) with CsCl-type B2-type structures. The calculations have 
been carried out the density functional theory and density functional perturbation theory in combination 
with the quasiharmonic approximation. The phonon-dispersion curves and phonon total and partial den- 
sity of states have been investigated. Our results show that the contribution of RE atoms is dominant in 
phonon frequency, and this character agrees with the previous discussion by using atomistic simulations. 
The temperature dependence of various quantities such as the thermal expansions, bulk modulus, and the 
heat capacity are obtained. The electronic contributions to the specific heat are discussed, and found to be 
important for the calculated MgRE intermetalHcs. 
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1 Introduction 

Magnesium (Mg) alloys have been attracting much attention and especially attractive for the applications of aero- 
nautical and automotive industry because of their strength to weight ratio 111. However, the traditional Mg alloys often 
show low strength and creep resistance at high temperature and this feature is a very serious problem of limiting their 
applications. Recently, it has been reported that some rare-earth-Mg intermetallic compounds MgRE (where RE indicates 
a rare-earth element) with B2 structures (CsCl-type structure) have good creep and high temperature strength El El SI. 
So MgRE intermetalHcs are extremely attractive structural materials for applications in automobile parts and aerospace 
industries, while various studies have been undertaken of the magnetic properties, linear and nonlinear elasticity, stacking 
fault, and thermal properties for the B2-MgRE intermetalHcs iS O H [S [9l [IHl [lU [HI. Recently, Wu et al 1 13] have 
performed the atomistic simulations to study thermodynamical properties of MgY, MgDy, and MgPr by using the mod- 
ified analytic embedded atom method and their results are generally in agreement with the experimental data and other 
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theoretical results. However, since the ionic degrees of freedom are treated classically in their study, this simulation are 
not valid at temperatures comparable to or lower than the Debye temperature. In addition, the first principles calculation 
indicates that the Fermi energy occurs above a peak in the electronic density of states (DOS) and B2-MgRE intermetallics 
have the large electronic DOS at Fermi energy level |7 |. Therefore, thermal electronic excitations contribution of the 
thermal properties must be considered. The exact calculations of thermodynamic properties for intermetallic compounds 
of Mg with rare earth metals, especially both the phonon and electron contributions being considered, have received little 
attention. 

It is obviously crucial to take into account the thermodynamic properties such as thermal-expansion coefficients, 
specific heat at constant volume and constant pressure , and temperature-dependent bulk modules in the investigation of the 
properties of MgRE intermetallics as well as their applications. A simplified method for thermal expansion calculations 
is in the framework of density-functional theory (DFT) with a Debye-Griineisen based model lfT4l . The more accurate 
approach is first-principles DFT treatments within the quasiharmonic approximation (QHA) which provides a reasonable 
description of the thermodynamical properties of many bulk materials below the melting point (TS] [T6| [iTl [TSl [T9ll . 
Quasiharmonic approximation lets one take into account the anharmonicity of the potential at the first order: vibrational 
properties can be understood in terms of the excitation of the noninteracting phonon. The exact calculations of phonon 
frequencies in every point of the Brillouin zone can be achieved by the density functional perturbation theory (DFPT) |[2Qll . 
Nowadays, quasiharmonic approximation based on DFPT has been applied with great success to more and more complex 
materials such as alloys [NiAl and NiAls EH], perovskite [MgSiOs EJl], and hexaborides [LaB6 and CeB6 |23 1]. In this 
paper, we apply first principles calculations within quasiharmonic approximation to study the thermodynamic properties 
of light RE MgPr, heavy RE MgY, MgDy, and MgTb intermetallics with B2 structures. The Phonon-dispersion curves and 
phonon density of states have been discussed. Thermal expansions, temperature dependence of isothermal bulk modules, 
and heat capacities at constant volume and constant pressure are presented. In addition, we have considered the electronic 
contribution to the thermal properties, and our results show that the eff'ects of thermal electronic excitations for MgRE 
(RE=Y, Dy, Pr, Tb) are remarkable. 

2 Theory 

The equilibrium lattice volume V of a crystal with cubic symmetry at any temperature T, is obtained by minimizing 
the Helmholtz free energy, F, of a system. The free energy at temperature T and constant volume V is give by 

F{V, T) = Eo{V) + FM T) + Fei(K T), (1) 
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where Eq(V) is the ground state (T = OK) total energy of the crystal, /^vib(K T) is the vibrational free energy which is 
come from the phonon contribution. In the quasiharmonic approximation FyibCK T) is written in the form 

qA qA 

Here, the sum is over all three phonon branches A and over all wave vectors q in the first Brillouin zone, is the 
Boltzmann constant, H is the reduced Planck constant, and ojqAiV) is the frequency of the phonon with wave vector q and 
polarization A, evaluated at constant volume V. The first and second terms of Eq. ^ represent the zero-point (F^p) and 
thermal energies (F^^), respectively. For the thermal electronic contribution to free energy Fei = Eei -TSei- The electronic 
energy E^i due to thermal electronic excitations is given by 

n(s,V)f(s)sds- n(s,V)sds, (3) 
Jo 

where n(s, V), f(s), and sf represent the electronic density of state (DOS), the Fermi-Dirac distribution, and the Fermi 
energy, respectively. The electronic entropy S ei is formulated as 



5ei(V, T) = -kB n(s, V)[minf(s) + (1 - /(^))ln(l - m)]ds. 
Jo 



(4) 



Usually, it is assumed that the electronic contribution to total free energy can be negligible. 

To calculate free energy, one must be able to calculate frequencies all over the Brillouin zone, and this can be done 
exactly using DFPT calculations. Furthermore, this calculation must be performed at various values of lattice parameter 
a in the case of isotropic thermal expansion for B2-structure crystal. The equilibrium lattice constant at temperature T is 
obtained by minimizing F with respect to a. The coefficient of volume thermal expansion is given by 



1 ( dV(T) \ 



V\ dT I 

and the linear thermal expansion is described by 

a(T) - a(Tc) 

^(T) = , (6) 

a(Tc) 

where a(Tc) is equilibrium lattice constant a(T) = [V(T)]^^^ at Tc = 300K. 

Once phonon spectrum is obtained, we can easily compute the temperature dependence of the vibrational specific 
heat capacity Cy and the entropy S at constant volume, as described below 

qA 

and 

/ no)^,(v)\ no)^,(v) no)^,(v)^ 

(8) 



= -kn}^[ in (2 smh - coth 

qA 



respectively. In addition, the electronic heat capacity can at constant volume be obtained from 




and we denote total specific heat at constant volume is then Cy = + Cy . 

Cp, the specific heat at constant pressure, can then be computed by using relation 

= Cy + a^BVT, (10) 

where B(T) = -l/Vd^F/dV^ is the bulk modulus. 

3 Computational details 

The first-principles calculations within the density functional theory (DFT) were performed using the VASP program 
ElEilESl. We employed the projector augmented wave (PAW) method (271 [28l. The Perdew-Burke-Emzerhof (PBE) 
1291 [3Q1 exchange-correlation functional for the generalized-gradient-approximation(GGA) was used. The structures are 
relaxed without any symmetry constraints with a cutoff' energy of 600eV for all calculated intermetallics. The Brillouin 
zones of the unit cells are represented by Monkhorst-Pack special k-point scheme fSTl with 21 x21 x21 grid meshes. 
The radial cutoff's of the PAW potentials of Mg, Y, Dy, Pr and Tb were 1.41, 1.72, 1.58, 1.61, and 1.58 A, respectively. 
The 3^ electrons for Mg, the 4^, 4p, Ad and 5^ electrons for Y, the 4/ and 6^ electrons for Dy, Pr, and Tb were treated as 
valence and the remaining electrons were kept frozen. In order to deal with the possible convergence problems for metals, 
a smearing technique is employed using the Methfessel-Paxton scheme |32|, with a smearing with of 0.05eV. 

We have carried out supercell approach within the framework of force-constants method for the phonon calculations. 
Real- space force constants of the supercell are calculated in the density-functional perturbation theory (DFPT) imple- 
mented in the VASP code 1331 from forces on atoms with atomic finite displacements, and the phonon frequencies are 
calculated from the fore constants using the PHONOPY package (341 [35l . From the full phonon spectrum, the lattice 
vibration free energy are calculated. Since the chosen supercell size strongly inffuences on the thermal properties, we 
compare the vibrational free energies of 3 x 3 x 3 supercell with those of 4 x 4 x 4 supercell at 300K and lOOOK, and find 
that the energy ffuctuations between 3x3x3 and 4x4x4 supercells are less than 0.01%. So the adequate supercell size 
consisting of 3 x 3 x 3 unit cells is chosen to calculate thermal properties for all calculated intermetallics; however it has to 
be considered that Grabowski et al 1 36 1 reported that small error in phonon calculation may bring significant change in the 
calculated thermal properties, especially at high temperatures. The thermal electronic energies and entropies are evaluated 
using one-dimensional integrations from the self-consistent DFT calculations of electronic DOS using FD smearing as 
shown in Eqs. ([3| and Q. In order to get the temperature dependence of lattice parameters, we have calculated total 
free energy at temperature points with a step of IK from to lOOOK at 15 volume points. At each temperature point, the 
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equilibrium volume Vt and isothermal bulk moduli B(T) are obtained by minimizing free-energy with respect to V from 
fitting the integral form of the Vinet equation of state (EOS) 1 37 1 at = 0. While, Cp has be calculated by polynomial 
fittings for Cy and by numerical diff'erentiation for dV/dT to obtain a(T). 

4 Results and discussions 

4.1 Phonon-dispersion curves and phonon density of states 

The calculated phonon-dispersion curves of four intermetallic compounds, which were calculated using 3x3x3 
supercells, are displayed in Figure [T] The phonon properties of all calculated intermetallic s are computed within the 
generalized gradient approximation (GGA) in the B2-type structure, with space group symmetry Pm3m(221), in which 
the Mg atom is positioned at (0, 0, 0) and the RE (RE=Y, Dy, Pr, Tb) atoms at (0.5,0.5,0.5). These materials contain 
two atoms per primitive cubic unit cell. Due to the symmetry, the dispersion curves are showed along high-symmetry 
direction F - X - M - F - 7? of the Brillouin zone. These dispersion curves have common framework. There is a phonon 
band gap starting around 3 THz. The maximum value of acoustic modes for MgY is slightly greater than 3 THz and 
those for MgDy, MgTb, and MgPr are slightly lighter than this value, since the atomic mass of Y is smaller than those of 
Dy, Tb, and Pr. The phonon density of states (DOS) including the partial DOS (PDOS) and the total DOS (TDOS) are 
shown in Figure |2] Sampling a51 x51 x51 Monkhorst-Pack grid for phonon wave vectors q is found to be sufficient in 
order to get the mean relative error in each channel of phonon DOS. The flat regions of phonon-dispersion curves, which 
correspond to the peaks in the phonon PDOS, indicate localization of of the states, i.e., they behave like "atomic states" 
tT9i . For the four intermetallic s, we find that the DOS are mostly composed of Mg states above the phonon band gap 
since its atomic mass is lighter than those of the rare earth elements RE (RE=Y, Dy, Pr, Tb). Our results agree with the 
previous theoretical investigation in which the contribution of RE atoms is dominant in phonon frequency since they are 
heavier than Mg atoms 11131 . 

4.2 Bulk properties and thermal expansion 

The present results of the equilibrium lattice constants gq and isothermal bulk modulus Bq at T = OK for MgRE 
(RE=Y, Dy, Pr, Tb) intermetallic s together with the available experimental values 1381 and the previous calculated results 
fTl are shown in Table [T] Our calculated results for the equilibrium lattice constants are within 0.8% of the experimental 
values, and shows excellent agreements with the previous theoretical results. For the bulk modulus which are obtained 
from Vinet equation of state, the present results are within 1.9% errors from the previous calculated values. The isothermal 
bulk modulus as a function of temperature B(T) are shown in Figure[3] Among the four intermetallic compounds, through 
the temperature range, heavy RE MgDy and light RE MgPr have the highest and lowest bulk moduli, respectively, i.e, 
MgDy is the most incompressible and MgPr is the most compressible. This indicates that the mechanical properties of 
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MgRE intermetallics can be further improved by addition of heavy RE metals |[T3l . With increasing temperature, the bulk 
moduli of the four intermetallics decrease and the differences among the bulk moduli almost remain unchanged. This 
character demonstrates that rare-earth-magnesium MgRE intermetallics can keep their mechanical properties through 
wide temperature range and have good high temperature stability, and it agrees well with recent reports that the MgRE 
have good strength at high temperature in comparison with the traditional Mg alloys EKJISI. 

The temperature dependence of the linear thermal expansion e defined by Eq. ^ of the four MgRE intermetallics 
are shown in Figure |4] The linear expansions of the four compounds are found MgTb>MgPr>MgDy>MgY, however the 
differences among them are small. The coefficients of the volume thermal expansion (2 as a function of temperature are 
shown in Figure[5]With increasing temperature, the thermal expansion expansion grow rapidly up to ~ 400^, and the slops 
become smaller and nearly constant at high temperatures. Below ~ 200^, those of MgDy, MgTb, and MgPr are equivalent 
and are greater than that of MgY. The reason can be understand from the electron configurations of the corresponding 
elements. We have the electron structures with Mg(3s^), Y(4d^5s^), Dy(4f^6s^), Fr(4f6s^), and Th(4f6s^). All the 
elements have two s electrons, and the outermostshell of rare-earth element Y includes the d electrons while the other 
elements with / electrons. So the thermal expansion of MgY is smaller than those of MgDy, MgPr, and MgTb. At high 
temperatures, those of MgTb and MgPr tend to be equivalent. 

4.3 Specific heat 

Once the phonon spectrum over the entire Brillouin zone is available, the vibrational heat capacity at constant volume 
Cy^ can be calculated by Eq. (|7j), while the electronic contribution to heat capacity at constant volume Cy can be obtained 
from the electronic DOS by using Eq. ([9]). Then, the specific heat at constant pressure Cp can be computed by Eq. ([TO]). 
As a comparison, both Cy values, including electronic contribution and not, are plotted. We display the results in Figure|6] 
As temperature increases, Cy^ tends to the classical constant 6R, and Cy and Cp still increase. At low temperature about 
below ~ 200^, the discrepancy between Cp and Cy can be neglected. For thermal electronic contributions to specific 
heat, we find Cy is not negligible at high temperature though smaller than Cy^. This character can be understood from 
the electronic DOS at the Fermi energy level fisp). For MgRE, they have large electronic DOS near Fermi energy level. 
When we consider the electronic contribution, the theoretical Cp keeps positive slopes obviously. Cy is smaller than C^^^, 
but it is larger than the value Cp - Cy = a^BVT. Usually, the electronic heat capacity can expressed as Cy = yT, where 
y is known as the electronic constants. To be guide to eye, we should mention that the order of these quantities from 
lower to higher is ^(MgPr) < ^(MgDy) < ^(MgTb) < ^(MgY). As previous studied by Wu et al |7|, the Fermi energy 
level for MgRE occurs above a peak in the electronic DOS; the bonding states are full, and filling of the anti-bonding 
states is sensitive to deviations in the local structure that affect the Fermi energy, so the electronic excitations affecting 
the thermal properties is remarkable. From Figure |6] it is clear that the electronic contribution should be considered for 
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MgRE intermetallics. 

5 Conclusions 

In conclusion, thermal properties of intermetallics MgRE (RE=Y, Dy, Pr, Tb) with B2 structures, such as thermal 
expansions, bulk modules, and heat capacities at constant volume and constant pressure as a function of temperature 
are studied by using the density functional theory and density functional perturbation theory in combination with the 
quasiharmonic approximation. It is found that the contribution of RE atoms is dominant in phonon frequency since Mg 
atomic mass is lighter than those of the rare earth elements RE. Through the temperature range, heavy RE MgDy is 
the most incompressible and light RE MgPr is the most compressible. With increasing temperature, the bulk moduli of 
the four intermetallics decrease and the differences among the bulk moduli almost remain unchanged. In the thermal 
expansion, we find MgTb>MgPr>MgDy>MgY at hight temperature, and below ~ 200K those of MgDy, MgTb, and 
MgPr are equivalent and are greater than that of MgY. The electronic electronic contributions to the specific heat are 
discussed, and found to be important for the calculated MgRE intermetallics. 
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Table 1: The equilibrium lattice constants ao and bulk modulus 5o at T = OK for MgRE (RE=Y, Dy, Pr, Tb) in 
our calculation in comparison with the previous calculated results and the experiment. 
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^ReferenceUJ; 

^Reference |[38ll , experiment. 
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Figure 1 : (Color online) Phonon-dispersion curves of (a) MgY, (b) MgDy, (c) MgPr, and (d) MgTb. 
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Figure 2: Phonon total density of states (TDOS) and partial density of states (PDOS) of (a) MgY, (b) MgDy, 
(c) MgPr, and (d) MgTb. The PDOS indicates that the density of states are mostly composed of Mg states at 
high frequency and RE states at low frequency. 
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Figure 5: (Color online) The coefficients of volume thermal expansion a as a function of temperature. 
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Figure 6: Temperature dependence of heat capacity of (a) MgY, (b) MgDy, (c) MgPr, and (d) MgTb. Our results 
shows that the thermal electronic contributions to specific heat is remarkable at high temperature. 
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